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phonon component is the decisive factor in the elementary excitations. We analytically solve the 
half-space boundary value problem of the temperature jump at the boundary of the degenerate 
Bose gas in the presence of a Bose — Einstein condensate. 

Keywords: degenerate quantum Bose gas, collision integral, Bose — Einstein condensate, 
phonon component, temperature jump, Kapitsa resistance. 



PACS numbers: 05.20.Dd Kinetic theory, 05.30.Jp Boson systems, 05.60.-k Transport processes, 03.75.Nt 
Other Bose — Einstein condensation phenomena 



1. INTRODUCTION 

The behavior of quantum gases has aroused increased 
interest in recent years. In particular, this is related 
to the development of experimental procedures for 
producing and studying quantum gases at extremely low 
temperatures The bulk properties of quantum gases 
have been studied in the majority of papers Ql and 

At the same time, it is obviously important to take 
boundary effects on the properties of such systems into 
account. We mention a paper where the thermodynamic 
equilibrium properties of quantum gases in a half-space 
were considered [4|. 

Along with the equilibrium properties, the 
nonequilibrium properties of quantum gases bounded in 
space attract interest. In particular, such a phenomenon 
as the temperature jump [5] at the interface between 
a gas and a condensed (in particular, solid) body in 
the presence of a heat flux normal to the surface is 
important. Such a temperature jump is frequently called 
the Kapitsa temperature jump [6|. 

We note that up to now, the Kapitsa jump has been 
calculated in the regime where only phonon scattering at 
the interface between two media was taken into _account 
and phonon scattering in the bulk was neglected [7]. Here, 
we take namely the effect of scattering of elementary 
excitations of the Bose gas (phonons) in the bulk into 
account. 



We take the character of phonon scattering by the 
surface into account by introducing a phenomenological 
coefflcient of specular phonon scattering by the surface. 
This approach is thus additional to that proposed in 

In [8|, we considered the temperature jump problem in 
a quantum Fermi gas. We obtained an analytic solution 
for an arbitrary degree of gas degeneracy. In we 
considered a similar problem for a Bose gas. But the 
gas was assumed to be nondegenerate, i.e., we did not 
take the presence of a Bose — Einstein condensate into 
account. 

This paper is devoted to analyzing the temperature 
jump problem in a degenerate Bose gas. The presence of 
a Bose — Einstein condensate leads to a considerable 
modification of both the problem statement and the 
solution method. In this case, to describe kinetic 
processes near the surface, we use a kinetic equation with 
a model collision integral. We assume that the boundary 
conditions at the surface are specular-diffuse. 

2. DERIVATION OF THE KINETIC EQUATION 

To describe the gas behavior, we use a kinetic equation 
with a model collision integral analogous to that used to 
describe a classical gas. We take the quantum character 
of the Bose gas and the presence of the Bose — Einstein 
condensate into account. For a rarefied Bose gas, the 
evolution of the gas particle distribution function / can 
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be described by the kinetic equation 

df as, 



dt dp 



V/ -/[/], 



(1) 



where £ is the kinetic energy of gas particles, p is the 
gas particle momentum, and /[/] is the collision integral. 

In the case of the kinetic description of a degenerate 
Bose gas, we must take into account that the properties of 
the Bose — Einstein condensate can change as functions 
of the space and time coordinates, i.e., we must consider 
a two-liquid model (more precisely, a "two- fluid" model 
because we consider a gas rather than a liquid). We let 
Pc ~ Pc{^,t) and Uc = Uc(r,t) denote the density and 
velocity of the Bose — ^instein condensate. We then can 
write the expressions 



n, 



PcUaUck- 



for the densities j , Q , and 11^^ of the respective mass, 
energy, and momentum fluxes of the Bose — Einstein 
condensate (under the assumption that the chemical 
potential is zero). The conservation laws for the number 
of particles, energy, and momentum require that the 
relations 



vn 



e{p)I[J]dnB. 



be satisfied in the stationary case. Here, 

(2s + l)(fip 



s is the molecule spin, £(p) is the energy, h is the 
Planck constant and /(/) is the collision integral in Eq. 
(l). 

In what follows, we are interested in the case of 
motion with small velocities (compared with thermal 
velocities). We note that for the Bose — Einstein 
condensate, the quantities Q and 11^^ depend on the 
velocity nonlinearly (they are proportional to the third 
and second powers of the velocity). Therefore, in the 



approximation linear in the velocity Uc , the energy and 
momentum conservation laws can be written as 



J e{p)I[f]dnB - 0, 



/ 



pi[f]dnB - 0. 



According to the Bogolyubov theory, the relation for 
the excitation energy £(p) holds for a weakly interacting 
Bose gas 



e{p) 



2 2 , 
UqP + 



P 

2m 



1/2 



(2) 



where 



Mo 



1/2 



a is the scattering length for gas molecules, n is the 
concentration, m is the mass, and p is the momentum 
of elementary excitations. The parameter a characterizes 
the interaction force of gas molecules and can be assumed 
to be small for a weakly interacting gas. 

In our previous paper ll|, we considered the case 
where the relation 

, kT 
m 

holds for sufficiently small a , where k is the Boltzmann 
constant and T is the temperature. In that case, the 
first term in the brackets in (2) can be neglected. The 
expression for the energy £ {p) takes the same form as in 
the case of noninteracting molecules: 

^2 



p 

2m' 



We now consider the case where the phonon component 
dominates in (2), i.e., where 



In this case, we obtain 

£(p) = uoIpI = uop. 
according to relation (2). Consequently, 

5£(p) P 
dp p 

When considering kinetic equation (1), by gas 
molecules, we must understand the elementary 
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excitations of the Bose gas with energy spectrum 
(2). The character of the elementary excitations is 
manifested in the properties of the colhsion integral. 
As a collision integral in Eq. (1), we take its r- 
approximation. Then the character of the elementary 
excitations is manifested in the dependence of the 
collision rate on the excitation momentum [ol, [lo| . 

ni-Qi 



^0=-— = Kp)(/B-/)• 



5/ ^ pdf 
at p or 

Here, / is the distribution function 



(3) 



J^(p) = vo\p- Pop 

is the dependence of the collision rate on the excitation 
momentum, and 7 is a constant. In the case where 
the phonon component dominates in the elementary 
excitations, 7^8 [5]. We have po = rnvg , where vq 
is the velocity of the normal component of the Bose gas, 
fg is the equilibrium function of the Bose — Einstein 
distribution 

-1 



cxp 



- 1 



and is a model parameter having the meaning of the 
inverse mean free path I . vq ^ l/l . 

The parameters in /g, namely, and u, , can 
be determined from the requirement that the energy and 
momentum conservation laws 



'^(p)p f-fh d^P = 0, 



y Kp)e(p)[/-/^]rf^p = 



be applicable. The conservation law for the number of 
particles is inapplicable here because of the transition of 
a fraction of particles to the Bose — Einstein condensate. 

We now assume that the gas velocity is much less than 
the mean thermal velocity and the typical temperature 
variations along the mean free path / are small compared 
with the gas temperature. Under these assumptions, the 
problem can be linearized. 

We seek the distribution function in the form 



where 



/ = /^(C)+¥.(i,r,C)5(C), 



exp C — 1 ' 



<y9 is a new unknown function, Tg is the surface 

d 

temperature, g(v) = —-^ — Eg = C, and we 



introduce the notation 



and 



mqP 



mqIP - p., I 



We then have 



9iC) 



exp(£*) — 1 ' 



exp C — 1 ' 



exp C 



(exp C- 1)2 ■ 

We linearize the local Bose — Einstein distribution 
(Bosean) fg , passing to dimensionless quantities. We 
note that 

Ts wo|p-p*| Ts uq 



^(p-p*)2 



or 



/ 5T^\f CC 



C 



CC, 5T^ 



Therefore 



r _ CC, ^^Tf, 

0£* — £* — Es — ^ O 'Tj^j 



because — C. 
Consequently, 



We note that the quantity J^(p) in Eq. (3) can 
be replaced with vqP'^ iii tti6 approximation under 
consideration. 

We introduce the dimensionless time 

r = Vol t 

\ Uq / 



and coordinate 



fkTsy 

ri = t'o r- 
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It is now clear that Eq. (3) (in the dimensionless 

variables) becomes 



dip C dip 



(4) 



From the energy and momentum conservation laws, we 
find 



in what follows. Such a function satisfies the equation 

dh 



i OO 

K{i,,C;tx',C')h{x,ii',C')g{C')dii'dC', (6) 

-1 



= / C^+^ip{T,ri,C)g{C)d^C, where K{ii,C; n' ,C') is the kernel of Eq. (6), 

A-KQ^+i J 



7+3 



where 



OO 

g^+n = j C^+"g{C) dC, n = 0, 1, 2, • • • . 



We represent Eq. (4) in the form that is standard in 
transport theory: 



d(p C 
dr C dvi 



+ CM^,ri,C) = 



= k{C, C'Mt, ruC')g{C')d'C, (5) 

where k{C,C') is the kernel of Eq. (5), 

Equation (5) can be represented in the equivalent form 



J C'^C'ipg{C')d^C'+ 



+ 



3. PROBLEM STATEMENT 

In the problem under consideration, a degenerate Bose 

gas occupies the half-space a; > above a planar surface 
where the heat exchange between the condensed phase 
and the gas occurs. Therefore, the function (p can be 
regarded as 

(P{t, ri,C) = h{x,iJ.,C) 



K{|^,C■,^J,',C') = ^ ^ + . 

57+3 .97+4 

The problem is to find the value of the relative 
temperature jump 

AT 

^ s 

where 



AT = Ts- T, 



as a function of Qx , which is the projection of the heat 
flux on the x axis. Taking linearity of the problem into 
account, we can write 

St = RQx- 

The dimensionless coefficient R of the temperature 
jump is called the Kapitsa resistance. 

It is obvious that Eq. (6) has the particular solutions 

hi{x,^,C)^fi and h2{x, fi,C) = C, 

and the Chapman — Enskog distribution function is 

has {x,iJ,,C) = B+n- EtC, 

where the quantity is proportional to the heat flux 

Qx • 

Assuming that the reflection of the elementary 
excitations from the wall is specular-diffuse, we now 
formulate the boundary conditions 

h{0,fx,C) =qh{0,-i2,C), 0</x<l, (7) 
h{x,ii,C) = 

= B+fi- etC + o{l), x ^ +00, -1< IJ.<0, (8) 

where q is the specular reflection coefficient. 

The problem is to solve Eq. (6) with boundary 
conditions (7) and (8). Finding the value of the 
temperature jump et is of special interest. 
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4. REDUCTION TO THE INTEGRAL EQUATION 

We continue the function h to the half-space a; < 
symmetrically: 

h{x,fjL,C) = h{—x,—iJ,C), X < 0. 

For X < , we then have 

has {x,fi,C) = B'fJ.- StC, 

with B+ = -B- . 

We now separate the Chapman — Enskog distribution 
from the function h assuming that 

h{x, fx, C) = B^^ — etC + hc{x, fx, C). 

For the function hc{x,fj,,C), we formulate the 
boundary conditions for the lower and upper half-spaces: 

/le(+0,M,C) = 

= -(1 + q)B+n + (1 - q)eTC + qhd+O, -/l, C), 
where < fi < 1 , 

/ie(-0,M, C) = 

= -(1 + q)B-ix + (1 - q)eTC + qK{-0, -ijl, C), 
where — 1 < /x < , and 

/ic(+oo, ^, C) = 0, hc{—c>o,fi,C) = 0. 

We include these boundary conditions in the kinetic 
equation. We obtain the equation 

+(1 - q)sTC - (1 - g)/ie(TO, C)}S{x). (9) 
Here, i5(a;) is the Dirac delta function, and 



1 oo 



+ 



j j ii"'C^+^h{x,fj,,C)g{C)dfj,dC, m = 0,l. (10) 

-1 



Equation (9) actually combines two equations. The 
point is that the term ft,c(— 0,/i, C) corresponds to 
positive fi : < jj < 1 and the term /ic(+0, /k, C) 
corresponds to negative : — 1 < < . 

We seek the solution of Eqs. (9) in the form of Fourier 
integrals: 

oc 

hc{x,fx,C) = ^ J e^'="$(fc,M,C)dfc, 



OO 



OO 

Wo{x) = ^ J e"'''Eo{k)dk, 

— OO 

OO 

Wi{x) = ^ I e''"'Ei{k)dk. 
2-K J 



Solving Eq. (10) with x > and < , assuming 
that the right-hand side of this equation is known, and 
assuming that the boundary conditions are satisfied far 
from the wall, we obtain 

h+{x,fi,C) = 



+ 00 

= -^exp(-^C^) J exp(^CT)w(t,M,C)rfi, 



where 

After simple calculations, we obtain 

+ 00 



h+{x,ii,C) 



f e'''''E{k,^,C)dk 



2n 



CT + ikjj, 



where 



E{k,ii,C) = ^r^Eoik) + -^Ei{k). 



257+4 

We can similarly show that 



2^7+3 



C7 /• E{k, II, C)dk 



Ct + ikn 
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From the last two expressions, taking the evenness of 
the function E{k,fi,C) with respect to the variable k 
into account, we obtain 

/i±(0,/x,C) = 



two characteristic equations. Without writing them, we 
combine them into one vector characteristic equation 

A{k)E{k) = -2(1 + q)B+Ti{k) + 2(1 - q)eTT2{k)- 



+ 00 



CT r {C^ -ikn)E{k,tx,C)dk 
2tt J C2^TfcV 



C^T ^rE{k,fi,C)dk 



IT J C27 + A;>2 • 



(11) 



Hence, we see that two Eqs. (9) can be combined into 
one: 

fi dh 



dx 



hc{x,fj,, C) = 



|M|{-(l + g)fi+|M|- 



+ (1 - q)eTC - (1 - q)htiO,fi,C)}S{x). (12) 

We pass to the Fourier integrals in Eq. (12) and obtain 
the equation 

{C^ +iktxMk,f,,C) = Eo{k) + —JtE^{k)+ 

^9l+i ^57+3 



-^-^ j J{k,ki)E{ki)dki. 





(15) 



Here 



1 oo 



Tm,i{k) 







li'C"'g{C)dfj,dC 

C27 + ^2^2 ' 



A(A;) is the dispersion matrix 



m 



1 —Ts^+Afiik) -^^T2T,+3,2(fc) 

07+4 ^7+3 

ik 3 

727+4,2(fc) 1 13^+3,2 (fc) 

57+4 57+3 



E{k) is an unknown column vector 

m 



Eo{k) 
Ei{k) 



Ti{k) and T2{k) are the column vectors of constant 
terms 



Ti(fc) = 



Ts-y + sAk) 

-ifcT2^+3,4(fc) 



{l+q)B+\n\ + {l- q)sTC - C^-^ X 



CEoiki) ^ 3,iEiiki] 



2^7+4 ^9"/+3 



dki 



C27 + fc2^2 



(13) 



Ej{k) 



1 oo 



+ J J fx^^k,fx,C)C''+''giC)dfxdC, i = 0,l. (14) 

-1 



Mk) = 

-ZfeT2^+4,3(fc) 

and J{k,ki) is the matrix function 



We consider equality (10). Rewriting it in terms of the •^(^> ^i) 
Fourier integrals, we obtain 



J47+4,l(fc, /Ci) ^^^37+3,3(^,^1) 

.97+4 57+3 

-^^-^37+4,3(^,^1) —^^"47+3,3(^,^1) 
57+4 57+3 



which we call the Neumann matrix. 
The integrals 



1 00 



Jm,i{k,ki) 



H^C'^g{C)diJLdC 







(C27 + A;2/i2)(C27 + fc2^2) 



We solve Eq. (13) for ^{k,fi,C) and substitute it 
in (14) for j = 0, 1 . We obtain a system consisting of are introduced in the Neumann matrix. 
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5. METHOD OF SUCCESSIVE 
APPROXIMATIONS 

We seek the solution of Eq. (15) in the form 

1 + g 



£T = 



1-5 



and 



£o + ei(l -g) +£2(1 -g)' 



E{k) = 2{l + q)x 



(16) 



E^°^ (k) + £(1) {k){l-q) + ^(2) (fc) (1 _ g)2 + . . 1 (17) 



We substitute these equahties in the characteristic 
equation. We obtain the countable system of equations 

Aik)E''°\k) = -B+Ti{k) + eonik), (18) 



A{k)E^^^{k) = £iT2(fc) - J{k,ki)Eo{ki)dki, (19) 



k%+4,2{k)T^+3,4{k) + T2^+3,2{k)T2j+4,2{k) 

The matrix inverse to the dispersion matrix is 

Dik) 



A-\k) 



A(fc)' 



, where 



Dik) 



Tn+3AW J2n+3,2(K) 

5ra+3 9n+3 



ik 



<-2n+i,2 



(k) 



-Tn+i,2{k) 



9n+4 571+4 

We consider the construction of series (17). We set 



E(°\k) = 



Ei'Hk) 



We consider the construction of series (17). We set 



A{k)E^^\k) = 

OO 

{k)-l J J{k,ki)Ei{ki)dku--- 





(20) 



We transform the elements on the principal diagonal 
of the dispersion matrix: 

1 k'^ 

1 Ts^+ifiik) = T^+4,,2{k), 

9i+i 9j+4 



9j+3 



'37+4 



9-i+z 



The dispersion matrix function now becomes 



-T^+A,2{k) 



3ik 



'27+3 



,2{k) 



9^+4 57+3 

--l27+4,2(fc) 7-i7+3,4(K) 



57+4 



9l+3 



We call the determinant of the dispersion matrix the 
dispersion function: 



X{z) = det A(z) = k^ujik), 



where 



uj{k) 



ff7+3ff7+4 



-B+T3^+3,2(fc) + eoT37+4,l(fc) 
ik{B+T2^+3A{k) + £o727+4,3(fc)) 



(21) 



Substituting expression (21) in Eq. (19), we find 
E^'^\k) . Substituting E^^\k) in Eq. (20), we then find 
E^'^\k) . Cont inuing this process without bound, we 
construct all terms of series (17). 

We now show how to construct the terms of series (16). 
For this, we write vector equality (21) in the form of two 
scalar equalities: 



E, 



(0) 



k'^uj{k) 



9-1+3 



'7+3,3 



(fc)x 



S+T3^+3,2(fc)-£oT3^+4,lW 



+ 



+ T2^+3,2{k) S+r2^+3,4(fc)+ £0^27+4,3 (A:) ) (22) 

07+3 L J 



(0) 



fc2w(A;) 



|ifcr2^+4,2(A;) B+T3^+3,2(fc) - £023^+4,1 (A:) 



+ 



+ik''T^+4,2{k) [B + T2^+3,4(fc) + £0^2^+4,3 (fc)] }• (23) 
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It can be seen from expression (22) that the function 
E^°\k) has no singularities at zero, while the function 
E[^\k) according to (23) has a second-order pole at A: = 
. Eliminating the singularity at zero, we obtain 



'37+4,1 



(0) 



3^7+4 



The quantity S+ is proportional to the heat flux: 

Q = 



/(x,p)^e(p)dnB. 



dp 

We transform this expression as 
(2s + l)u§ 



Q 



{2-Khf 

Integrating in this expression over the dimensionless 
momentum, we obtain 

(2s + l)(fcT,)4 



Q 



h{x, C)g{C)Cd'^C. 



We replace here the function h{x, C) with its 
Chapman — Enskog expansion has{x,C) . As a result, 
we have 



{2s + l){kTs) 



-10 



for the X component of the heat flux. 
We hence have 

67r2 



B+ 



53 (2s + l)(A:T,)4- 
The quantity eo is therefore 



So = 



(24) 



.93.97+4 (2S+I)(fcr,)4 

Returning to the formula for the temperature jump 
AT = RQ^, 

we find from expression (24) that the Kapitsa resistance 
in the zeroth approximation is 



R = 



Hh 



g^g^+i 1-9 (2s + l)fc4r3- 
We rewrite this equation in the form 



i? = C(7,9) 



{2s + l)k^Tf ' 



C(T,q) 14 




Fig. 1. Dependence of the temperature jump coefficient on 

the parameter 7 : the specular reflection coefficient is 
q = 0.3 for curve 1,9 = 0.5 for curve 2 , and 9 = 0.8 for 
curve 3 . 



C(r,q) 




Fig. 2. Dependence of the temperature jump coefficient on 
the specular reflection coefficient q: the collision parameter 
is 7 = 3 for curve 1 , 7 = 4 for curve 2 , and 7 = 5 for 
curve 3 . 



where 



, (25) 

ff357+4 1-9 

is the (dimensionless) coefficient of the temperature 
jump. 

The plus sign in formula (25) indicates that the 
wall temperature is higher that the phonon component 
temperature. 
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The graphs of the behavior of the temperature jump 
coefficient as a function of the parameter 7 and the 
specular reflection coefficient q are shown in Figs. 1 and 
2. 

It can be seen from these graphs that the quantity 
C{'y,q) decreases monotonically as the parameter 7 
increases. As the specular reflection coefficient q tends 
to unity, the quantity C{'y,q) increases without bound 
because the heat exchange between the wall and the gas 
adjoining it becomes impossible in this limit. 

6. CONCLUSIONS 

We have constructed a kinetic equation for a 
degenerate quantum Bose gas whose collision rate 
depends on the momentum of elementary excitations of 
the Bose gas. We considered the case where the phonon 
component is the key factor in the elementary excitations. 
The boundary conditions were assumed to be specular- 
diffuse. We solved the half-space boundary value problem 
of the temperature jump at the boundary of a degenerate 
gas in the presence of a Bose — Einstein condensate. We 
derived a formula for finding the temperature jump and 
calculating the Kapitsa resistance. 

We have developed a sufficiently general method 
for solving kinetic equations with specular-diffusive 
boundary conditions; this method was first proposed in 
where the problem of the skin effect was considered. 
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